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Abstract 

We use the observed SU (3) breaking in the mass spectrum of mesons contain- 
ing a single heavy quark to place restrictions on the light quark current masses. 
A crucial ingredient in this analysis is our recent first-principles calculation 
of the electromagnetic contribution to the isospin- violating mass splittings. 
We also pay special attention to the role of higher-order corrections in chi- 
ral perturbation theory. We find that large corrections are necessary for the 
heavy meson data to be consistent with m u = O.We use the observed SU (3) 
breaking in the mass spectrum of mesons containing a single heavy quark to 
place restrictions on the light quark current masses. A crucial ingredient in 
this analysis is our recent first-principles calculation of the electromagnetic 
contribution to the isospin-violating mass splittings. We also pay special at- 
tention to the role of higher-order corrections in chiral perturbation theory. 
We find that large corrections are necessary for the heavy meson data to be 
consistent with m u = 0. 
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I. INTRODUCTION 



The light quark current masses m u , m<f, and m s are among the fundamental parameters 
in the standard model of particle interactions, but an accurate and reliable determination 
of these parameters remains elusive. The reason is that these masses are small compared 
to the mass scale associated with confinement, A ~ 1 GeV, so that the light quarks are 
tightly bound inside hadrons and their mass cannot be measured directly. However, we 
can expand around the chiral limit m u ^,s — * to obtain information about current mass 
ratios using chiral perturbation theory |l]-(3[]. The idea is that in the chiral limit, there is a 
SU (3) l x SU (3)r chiral symmetry that is spontaneously broken to SU(3) l + r at the scale A. 
In this limit, the theory contains 8 Nambu-Goldstone bosons, which are identified with the 
light pseudoscalar meson octet (n, K, rj). The low-energy interactions of these states can be 
parameterized by an effective lagrangian with a few unknown parameters In addition, 
the chiral symmetry is explicitly broken by the light quark masses and by electromagnetism. 
This breaking can be treated perturbatively in the quark masses and the electromagnetic 
coupling, and selection rules associated with the chiral symmetry again tightly constrain the 
form of these perturbations. 

The main difficulty in this approach is that m s /A ~ 20-30%, so that higher-order effects 
can significantly change low-order results. This is the source of the interesting (and contro- 
versial) issue of whether higher-order corrections can be large enough to allow m u = |5|-[7| , 
thus solving the the strong CP problem. This will be a large part of the focus of this paper. 

In this paper, we apply the methods of heavy quark symmetry [|8| p~0|| and heavy hadron 
chiral perturbation theory |TT| to the determination of the light quark masses from the mass 



spectrum of mesons containing a single heavy quark. These states have quantum numbers 
Pq ~ QQi where Q = b,c and q — u,d, s. If we could ignore electromagnetic effects and 
higher-order corrections in the quark masses, we would obtain 

= m s -m ? P s -P _ f 20 for P = D (1 ^ 

~ m d -m u P d - P u ~ { 280 for P = B [ ' ' 

where rh = \{m u + rrid), etc. (We use the names of the heavy- meson states to denote their 
masses.) The enormous discrepancies between these numbers clearly cannot be accounted 
for by the higher-order quark mass corrections. The reason for these discrepancies is simply 
that the electromagnetic contribution to Pd — P u is numerically comparable to the contri- 
bution from the light quark masses. This underlines the importance of determining the 
electromagnetic contribution to the mass differences. 

We will carry out an improved analysis of R making use of the heavy-meson electromag- 
netic mass differences computed in ref. ||12|| . (For earlier related work, see ref. |13| .) Our 
analysis will also include important non-analytic corrections in the quark masses p^JT5[] 
and will pay careful attention to the role of higher-order corrections in chiral perturbation 
theory. We will briefly describe the elements of the calculation of the electromagnetic mass 
differences in the next section. The following section contains our analysis of R, and the 
final section gives our conclusions. 
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II. COMPUTATION OF ELECTROMAGNETIC MASS DIFFERENCES 



In ref . [0 , we computed the electromagnetic mass differences of heavy mesons in terms 
of measurable data using techniques similar to those used in the classic calculation of the 
7r + -7T° mass difference |16| . The basic idea in both calculations is to use dispersion-theoretic 



arguments together with the ultraviolet properties of QCD to relate the electromagnetic mass 
difference to measured strong-interaction matrix elements. We review the main features of 
our calculation here in order to make this paper more self-contained. 
We begin by writing 

(EM) ie>_ r d 4 k AT(p,k) 



where AT = T d — T u is a difference of the (traced) Compton amplitudes 

T q (p,k) = i j d A xe ikx (P q (p)\J»(V)J,{x)\P q (p)), (2.2) 

where J M is the electromagnetic current. It can be shown rigorously that this integral 
converges [[H],|r3 • 

We consider AT in the large- N limit, where N is the number of QCD colors. In this 
limit, we can express AT as a sum of tree graphs in a theory of (infinitely many) mesons with 
interactions that are at most polynomial in momenta . (This can be viewed as a rigorous 
version of the "resonance dominance" assumption made in refs. ||16|| .) The vertices of the 
meson graphs that determine AT are directly related to physical processes. For example, 
the graphs in fig. lb determine the heavy-meson form factors, while the last two graphs in 
fig. la are related to meson scattering amplitudes. Thus, the representation of AT in the 
meson theory can be viewed as a double dispersion relation in k and p. 

We then impose constraints on the couplings appearing in this sum by demanding that 
the matrix elements have the infrared behavior required by electromagnetic gauge invariance, 
and the ultraviolet behavior demanded by QCD. These constraints can be expressed as sum 
rules relating the couplings appearing in the sum. (Two of these sum rules are close analogs 
of the Weinberg sum rules that appear in the calculation of the 7r + -7r° mass difference [[H^].) 

Even after imposing the sum rules, the electromagnetic mass difference is given by an 
infinite (convergent) sum over intermediate states. We then derive an unsubtracted fixed-fc 2 
dispersion relation for AT that shows that the contributions of intermediate states with large 
mass are suppressed by inverse powers of their mass. It is therefore sensible to assume that 
the sum is dominated by the lowest-lying intermediate states. That is, we keep the minimum 
number of intermediate states to saturate the known asymptotic behavior of the form factors 
and other matrix elements that appear. This turns out to be the P and P* heavy-meson 
states (which are degenerate in the limit uiq — > oo) and the vector mesons p, u, p' and u>' 
that are responsible for "softening" the form factors. When the sum rules are imposed on 
this restricted set of intermediate states, there are sufficiently few free parameters that we 
are able to obtain numerical estimates for the electromagnetic mass differences. The results 
are 
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Here (/?') measure P*P r y (P*P* r y) couplings. Heavy-quark symmetry gives = (3' + 
0(1 /tuq), and experimental constraints from .D* decays give < — < 5 GeV -1 |20|j2l[| .PI 
Since the isospin splittings in the P and P* systems are equal in the heavy-quark limit, we 
cannot use them to get independent constraints on the quark masses. For the D system, 
we have chosen to work in terms of the "spin-independent" combination D + 3D* ||15|| , for 
which the the largest /9-dependent term in the electromagnetic mass difference cancels. We 
will use the D hyperfme splitting D* — D to give a determination of (3 below. 

There are a number of simultaneous approximations made in this computation: the 
heavy-quark expansion (including 1/mg terms), the large- N limit, and the truncation of 
intermediate states. The heavy-quark expansion appears to be under good control in our 
calculation, in the sense that the 1/mQ corrections are small compared to the leading terms. 
The large-iV limit and truncation of intermediate states are equivalent to approximations 
that can be made in the computation of the 7r + -7r° mass difference, which then works to 
30%. We conclude that it is reasonable to expect our estimates of the electromagnetic mass 
differences to be accurate to about 30%. 

We now complete our determination of the electromagnetic mass differences by extracting 
a value for (3 from the D hyperfme splittings. To order m q /mQ in chiral perturbation theory, 
one finds \Kl 
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Experimentally, (D* - £>*) - (D s — D) = 0.5 ± 2.5 MeV, and R > 20. Thus, the second 
term on the right-hand side is negligible, and the isospin hyperfme splitting is dominated 
by the electromagnetic contribution. Using eq. ( |2.5|) , we obtain 







■1.6 GeV" 1 



(2.7) 



x Our convention for (3 differs by a sign from that of ref. [20]. 
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for reasonable values of f3'. This is consistent with the results of refs. pOU2lf . For comparison 
the constituent quark model prediction is (3 ~ —3 GeV" 1 . It should be emphasized that our 
determination has errors 0(l/m c ), and there is a further 0(l/m c ) error in using this value 
for the B system. Using eq. ( |2.7|) , we obtain 



(EM) 



(D + - D°) + 3(D* + - D*°)\ v ""' ~ 2.5 MeV, (2.8) 
(B+ - B°) {EM) ~ 1.8 MeV. (2.9) 



The electromagnetic mass differences are not very sensitive to (3, so despite the uncertainties 
discussed above, we expect these values to be correct to about 30%. 

III. CONSTRAINTS ON R 

We now turn to the analysis of R defined in eq. ( |1 . 1| ) . We pay special attention to the 
role of higher-order corrections in chiral perturbation theory, since it is known that these can 
have an important effect on the determination of the light quark mass ratios |5j,|7]]. Keeping 
terms of order m q , m 3 / 2 , m g /mQ, m 2 , and the corresponding logs,0 the mass differences of 
the heavy mesons can be written in the form 

P d - Pu = (P d - Pu) iEM) + A (m d - m u ) + A^ + [AM + AM] m ^ m - 

m Q 

+ [A 2 (fx) + A 2 (/i)] m s (m d - m u ) + 0(m 2 u4 \nm s ) (3.1) 
P s -P = A m s + Ap) + [Ax(ai) + Ai(/i)] 

m Q 

+ [A 2 ^) + AM] < + [AM + AM)] ml + O(mf). (3.2) 

(For the D system, P = j(D+3D*).) The meaning of these terms is as follows: (Pd — -P U )^ EM ^ 
is the electromagnetic contribution discussed in the previous sections. Aom q is the term 
linear in quark masses that gives rise to eq. (|1 . 1|) . A^ 3//2 ) are the 0(m?J 2 ) non-analytic 



corrections, which can be expressed in terms of the light meson masses and couplings [TjJT 

AK? = -j£j- 2 (M 3 K0 - M| + + ^M 3 ), (3.3) 

Ai 3/2) = - T |^(Mi + iMj). (3.4) 

Here, / ~ 110 MeV is the K decay constant and g measures the strength of the PPn 
coupling [11]. The experimental limit on the D* width gives the bound g 2 < 0.5 [111. (To 



2 This corresponds to assuming that m s /A ~ A/rriQ. This probably overestimates the 0{m q /mQ) 
corrections for the B system, but this is harmless because our results are identical if these correc- 
tions are omitted. 

3 Ref. [Q] contains an error that was corrected in ref. [15|. 
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the order we are working, the value of g is the same in the D and B system by heavy quark 
flavor symmetry.) However, the analysis of refs. reveals that the magnitudes of g 

and (3 are correlated, so that g 2 ~ 0.15 for the value of f3 in eq. ( |2.7p . Their analysis also 
shows that g 2 > 0.1. We will conservatively assume 0.1 < g 2 < 0.3, which corresponds to 
-P < 3 GeV- 1 . 



The terms proportional to A% and A 2t 3 in eqs. (|3.1| , [3.2|) arise from analytic counterterms 
in the effective lagrangian of order m q /rriQ and m 2 , respectively. The terms proportional to 
A 12 3 ~ lnm s are the chiral logs that renormalize them. The counterterms and the chiral 
logs each depend on a renormalization scale [i in such a way that the terms proportional to 
Aj + Aj are independent of fi. The fact that the chiral logs are proportional to the same 
function of quark masses as the corresponding counterterms is a consequence of the simple 
structure of the SU(3) group theory for this system. 

Perhaps surprisingly, the chiral logs in A 2i 3 cannot be computed in terms of known 
quantities. The reason is the presence of higher-order terms in the effective heavy-meson 
lagrangian such as (in the notation of ref. p"T| ) 

8£ = T tr 
A 



tfttf (£tfl^ - _ £0 £t)l , (3.5) 



that are not constrained by present data. These terms give rise to non-analytic contributions 
such as 

A ' m -~!i^ ln f- 0(m ' lnm * ) - (3 - 6) 



(A similar fact was noted for baryons in ref. [22].) In principle, the log corrections are 
enhanced over the counterterm contributions by ~ In M\ / fi 2 for fi ~ 1 GeV, but in practice 
the logs are not significantly larger than unity. We will therefore treat the terms proportional 
to Aj + Aj as unknown corrections. 

We can solve for R from eqs. (|3.1| , |3.2j) to obtain 

R = (Ps-P)-^ /2) (K)-A 3 

(P d - P u ) - (P d - P u )(™) - A^l\K) ' 

where A^ 3 ^ 2 \K) are the terms in the 0(m^ 2 ) corrections which depending on the K 
masses (see eqs. (|3~3|j3^ )); A 3 = [A 3 (//) + A 3 (//)] m 2 parameterizes the unknown 0(m 2 ) 
and 0(m 2 \nm g ) deviations from the 0(m^ 2 ) mass relations for both the D and B systems. 
Since such relations are expected to work to 20-30%, it is reasonable to assume that A 3 
should not be much larger than 30%. Ref. || advocates a different measure of the chiral 
corrections. They demand that the second-order corrections to individual masses be less 
than 30%, but allow this to be the result of cancellations between larger corrections. Fol- 
lowing this criterion, we would allow A s to be 60%, since this can cancel against —30% 
corrections arising from the term proportional to A 2 + A 2 to give 30% corrections to P mass 
differences. 

We show R as a function of (Pd — P U Y EM ^ in figs. 2 and 3 for values of A3 corresponding 
to 0, 30%, and 60% of P s — P. We see that if the chiral corrections parameterized by A3 
are ~ 30%, and the computed electromagnetic mass differences have errors < 30%, then 
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the results for both the D and B systems can comfortably accomodate R = 44. This is the 
value of R predicted by lowest-order chiral perturbation theory for the light meson masses; 
it is also the value obtained by the O(mjj) analysis of ref. ||. 

We can use our results to address the question of whether m u = by using a relation 
between R and m u /md obtained from the light pseudoscalar masses that is valid to OimT) 
and that predicts R = 24 ± 2 for m u = 0. From figs. 2 and 3, we see that if we 
assume that A3 is a 30% correction, we require very large corrections in both the D and 
B electromagnetic mass differences in order to be consistent with m u = 0. If A 3 is a 60% 
correction, the data can accomodate m u = 0. 



IV. CONCLUSIONS 

We have derived constraints on the light quark mass ratio R from the spectrum of 
mesons containing a single heavy quark, using the QCD-based computation of the heavy- 
meson electromagnetic mass differences of ref. | 12| . Our results, summarized in figs. 2 and 



3, indicate that even when 30% uncertainties are assigned to both the electromagnetic mass 
differences and the unknown 0(m;j) chiral corrections, the value of R is bounded away from 
the value required by m u = 0. While we do not regard this as definitive proof that m u ^ 0, 
it is striking that the central values of higher-order analyses of both the light pseudoscalar 
mesons and the heavy mesons prefer m u ^ 0, and large chiral corrections must be invoked 
in both cases to be consistent with m u = 0. 
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FIGURES 



FIG. 1. Diagrams contributing to T in the large- N limit of QCD. The sum over n is over heavy 
mesons with quantum number Qq, while the sums over r and s are over light vector mesons with 
quantum numbers qq. 

FIG. 2. The quark mass ratio R as a function of the spin- independent electromagnetic mass 
difference in the D system. Our prediction for the electromagnetic mass difference is shown by the 
vertical dashed line. The bands correspond to the range 0.1 < g 2 < 0.3. The chiral corrections 
parameterized by ^3 are assumed to be 0% in the upper band, 30% in the middle band, and 60% 
in the lower band. The value R = 44 is the value obtained from a leading order analysis of the 
light meson masses, while the value R = 24 is the value required by m u = 0. 

FIG. 3. The quark mass ratio R as a function of the B electromagnetic mass difference. Our 
prediction for the electromagnetic mass difference is shown by the vertical dashed line. The bands 
correspond to the range 0.1 < g 2 < 0.3. The chiral corrections parameterized by ^3 are assumed 
to be 0% in the upper band, 30% in the middle band, and 60% in the lower band. 
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